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ABSTRACT 


A typical dead reckoning system utilizing a vehicle attitude 
and heading reference package is studied. The distance information 
is resolved into coordinate axes which yield vehicle position. System 
error models are developed which provide equations for position error 
versus time and distance. Constant velocity operation is assumed in 
the analysis however the error equations may be developed for any 
velocity profile where vehicle velocity is expressed as a function of 
time. A linear, single degree of freedom vertical loop, and a free 
directional gyro reference are the principal items studied in the error 
analysis. Component errors and vehicle disturbances are described 
and their effect presented in a mathematical relationship. 

An error analysis, where typical vehicle parameters are 
assumed, is shown, with results presented for both rotational errors 
of the coordinate system and for resultant position errors. 
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SECTION 1.0 


INTRODUCTION 


This report is a summary of the study effort during Task Order 
N-47 The Scope of Work was to include an analysis of a dead reck- 
oning navigation system for a lunar roving vehicle. To provide the 
required information a series of items are studied, relating to the use 
of a stabilized platform as a basic heading and attitude reference de- 
vice. The principal part of the report is a description of the system 
studied and an analysis of the errors obtained from the use of this 
system. Results show a cross track error of 300 meters and an alti- 
tude error of 70 meters for a mission time of six hours and for typical 
velocities of 5 to 10 kilometers per hour. 

The dead reckoning equipment is to be utilized by the vehicle driver 
during portions of lunar travel between position fixes. The guidance 
problem is not considered here inasmuch as guidance of the vehicle 
will be achieved primarily by visual observations for both the manned 
and the unmanned missions. The dead reckoning system is to be utilized 
intermittently by the driver to indicate his postion on a lunar chart. 

Tlis relation is used to provide an indication of his position in relation 
to the terrain, which is the primary pilotage data source. 



SECTION 2. 0 


SYSTEM REPRESENTATION 


A dead reckoning navigation system for a land vehicle utilizes 
distance and direction to compute present position in relation to a known 
starting point. The distance measurement is resolved into components 
along some selected set of axes. Since the driving of a vehicle along an 
arbitrary path over rough terrain involves frequent changes in heading, 
as well as pitch angle (as referred to the isogalic plane), it is necessary 
that distances along the coordinate axes be the result of an accumulation 
of incremental distances. 

The reference coordinate system will be a selenographic fixed 
location system as shown in Figure 1. The dead reckoning system will 
utilize starting coordinates consisting of a latitude and longitude reading. 
Distance traversed is to be related back to this starting point. In Figure 
1 the XYZ coordinates represent the inertial reference frame. The 
x tYt z t system represents the selenographic system in •which the vehicle 
navigates. The xtyt^t system is rotating about the Z axis of the inertial 
system with angular velocity SL . 

Now if we consider the vehicle moving on the surface of a 
sphere of radius R, which is assumed to remain constant, then the 
vehicle is located at any time by coordinates bl , rj where OC represents 
the latitude change and jO, the longitude change from the reference 
point. The mechanism for obtaining OL and yQis contained within the 
vehicle. CC and $ are of course, equal to xt/R and yt/R where: 

x^ is the distance along an arc of the great circle path in the 
direction the vehicle is traveling. 

y^. is the distance along the y direction, or the arc length of 
a direction transverse to the path of the vehicle. 

R is the radius of the body. 

The orthogonal coordiante system x , y^, z^ is constrained to the 
surface of the body hence, everywhere on tne body the z axis is directed 
along the radius vector. The system measures the rotation of the x , 
y , z t system from a inertial set at the starting point. Hence the object 
ox the dead reckoning system studied herein is to yeld the x , y , and 
z^ distances along the axes. Since these resolved distances are the 
result of an accumulation of incremental distances along the path they 
can be expressed mathematically as: 
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FIGURE 1 

REFERENCE COORDINATE SYSTEM 



where 


x t 


yt 

z t 


J( cos^)(cos p) ds 
J'(sin^JS) (cos p) ds 
^(sin p) ds 


* 


heading of vehicle with respect to north 


p = pitch of the vehicle relative to an isogalic plane 


ds = incremental path length 


One method of mechanizing these equations is shown in block diagram 
in Figure 2. Such a mechanization would utilize vehicle velocity as an 
input and perform the integration with respect to time. The outputs, 
when converted to angles and 0 can be used to provide the latitude 
and longitude coordinates. 


2. 1 MEASUREMENT OF ANGLES 

Assume that a coordinate system is established as shown in 
Figure 3. In this diagram x v is directed along the longitudinal axis of 
the vehicle, y v is along the traverse axis, and z v toward the vehicle's 
zenith. 


The relationship of the vehicle system to the selenographic, t, 
system is given by the Euler angles which relate the transformation 
between coordinate systems. This is shown in Figure 4. 

The vehicle relationship to the selenographic system is given 
by an azimuth angle, a, a roll angle, r, and a pitch angle, p. 

Figure 5 represents a line diagram of a gimbal arrangement 
for measurement of vehicle pitch angle with respect to the vehicle co- 
ordinate system. The vehicle will pitch and roll in its own coordinate 
system x v , y v , z v . To obtain the desired pitch angle with respect to 
the reference system x t , y t , z t the readouts of the synchros are fed 
through an appropriate transformation network. 

The measurement of the azimuth angle requires an additional 
gimbal. Azimuth is the angle between the x^. and the x & axis. A 
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DEAD RECKONING COMPUTER MECHANIZATION 











FIGURE 5 

TYPICAL GIMBAL SYSTEM 
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directional gyro or platform will normally be aligned toward north and 
measurements of the vehicle direction are made with respect to this 
reference. The directional gyro may be a separate instrument, or may 
be incorporated in a three axis platform. The required measurements 
are then the distance, the pitch angle, and the heading angle. If these 
quantities could be determined accurately then the mechanization of 
Figure 2 would yield the correct location of the vehicle. Errors in the 
measured quantities will cause the computed position to be in error, the 
total system error being the r. s. s. of the errors in measurement due 
to all the errors sources, each taken individually. The classification 
of errors to be considered is given in the next section. 
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SECTION 3. 0 


CLASSIFICATION OF ERRORS 


The three error source quantities d r , dP, and d^ may each 
be thought of as consisting of several factors which may be categorized 
into groups. Each error source is then to be classified as one of the 
following types. 

Systematic Error. An error that is a function of a system parameter 
susceptible to calibration and removal either at the initial alignment or 
at a later date. 

Noise Error . An error that occurs in a random fashion or the source 
of which is of a random nature. 

Bias Error. An error which is non-varying over short periods of time 
but which cannot be compensated for. 

The classification of the errors and the listing of errors sources 
for a navigation system is given in Table 1. 

Distance errors are not considered in this report inasmuch as 
information required to perform such an analysis is not available. For 
purposes of completness the estimate for accuracy of distance is taken 
as one percent of the distance traveled. Odometer systems in use by 
the U. S. Army are reported to have a system accuracy of this magni- 
tude (Reference 1). 
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TABLE 1 

CLASSIFICATION OF DEAD RECKONING 
SYSTEM ERRORS 


Type of Error 


Source of 
Errbr 

Systematic 

Bias 

Noise 

Distance 

Wheel diameter to meter 
Ratio 

Wheel distortion 

Wheel Slippage 

Initial 

Position 


Gravity 

Readout Operator 
Interpolation Tables 
Vertical Sensor 

Pitch Angle 

Velocity Error 
Moon Rate Error 
Coriolis 
Acceleration 
Vehicle Turn Rate 

Gravity 

Servo Dead Zone 

Vehicle Dynamics 
Readout Error 
Gyro Drift 
Vertical Sensor 
Linearity 

Azimuth 

Angle 


Initial Alignment 
Gyro drift 

Gyro drift 

Final 

Position 



Map Inaccuracies 
Display System 
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SECTION 4. 0 


SELECTION OF A SYSTEM 


The measurement of the pitch and roll angle of the vehicle is 
made with respect to a gyro held vertical, and a directional gyro aligned 
toward north. The reference direction for the vertical is the direction 
of a plum bob. 

A vertical directional reference is necessary for two reasons; 

(1) to provide an axis from which vehicle pitch may be measured, and 

(2) to maintain the spin axis of the directional gyro in a horizontal plane. 
Gyros are inertial instruments, therefore due to the rotation of the 
moon, and due to motion of the vehicle over the lunar surface, the 
inertial axes will not maintain their correct relationship with the seleno- 
graphic coordinate system. To maintain the directions over useable 
periods of time requires compensation methods. In the case of the 
vertical, the gyro spin axis can be made to maintain its direction along 
the axis by use of a vertical loop. The direction of the gravity vector 
is determined by a vertical sensing element, and the error measured 

is used to torque the vertical gyro. 

The north direction cannot be slaved to an external reference 
due to the slow spin rate of the moon and the absence of a magnetic 
field, and a north reference must be compensated for in an open loop 
manner. 


In order to study the accuracy of the system it is necessary to 
obtain equations which give the performance of both the vertical loop 
and the directional gyro. 

4. 1 OPERATION OF THE VERTICAL LOOP 

The vertical sensing element in an erection loop is an accelero- 
meter or a pendulum. Such a device of course is sensitive to vehicle 
accelerations as well as the acceleration due to gravity. If the pendulum 
on the stabilized element is tilted by some angle, 0, then the accelera- 
tion sensed by the dovice can be expressed as 

sV(s) - g0(s) 
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where 


s - Laplace Operator 

V - Vehicle Velocity 

g - Acceleration due to gravity 

The total angle through which the axis must be moved is equal 
to the tilt angle plus the angle due to rotation about the body of radius 
R. Since the total sensed acceleration should be made proportional to 
the rate of change of the correction angle, the following expression ma’ 
be written (see Ref. 2). 


K(sV(s)-g0(s) = 0 (s) + 

R 

0(s) (s+ K g ) = (Ks - 1 / Rs)V(s) 

0(s) _ Ks - 1/R 

V(s) £ + Kg 


Expressed in block diagram form 



FIGURE 6 

SIMPLIFIED LOOP DIAGRAM 


13 




Note that the erection loop has two inputs, one consisting of 
the vehicle acceleration sV(s), and the other the turn rate around the 


body 


V(s) 

Rs 


the first of these terms, sV(s), may be considered as a 


noise source. The bandwidth of the servo loop is to be adjusted so 
that this term, vehicle acceleration, is to be effectively filtered. This 
requires a reasonably low value of gain K. At the same time, due to 


input 


V(s) 

Rs 


a large value of gain is desired in order to reduce the 


tilt angle due to vehicle velocity. This tilt angle error is given by the 

0 { s \ 

final value theorem of the transfer function — i — — . 

V( s ) 


Let V(s) = — 
s 


1 t-*«> 

lim i \ Vo 

= sG s = 

s-»0 KgR 

where 


0 

« 

11 

steady state tilt error 

v 0 = 

velocity of vehicle 

K 

gain in the erection loop 

R 

radius of the body 


The adjustment of the gain K is then a compromise between 
filtering of vehicle induced perturbations and reduction of the velocity 
error. 


The velocity of the lunar vehicle is established between limits. 
A velocity of 20 ft/ sec is a possible upper limit, while much lower 
velocity levels will generally be the case. 
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Vehicle induced perturbations are those generated by accelera- 
tion of the vehicle along the direction of travel and transverse to the 
direction of travel, as caused by vehicle turn rates. It is difficult to 
establish values and times associated with this input, and more know- 
ledge of the vehicle and the terrain is necessary before this item can 
be determined. 

4. 2 SCHULER TUNED VERTICAL LOOP 


By proper choice of system gain, it is possible to reduce the 
velocity lag to a minimum, and make the system not susceptible to 
acceleration errors. This is obtained by the addition of a second inte- 
gration to the loop and the choice of gain K equal to 1/R. When these 
two conditions are met the vertical loop becomes Schuler tuned and 
the vertical direction, at least theoretically, is maintained independent 
of vehicle motion. 

The problems associated with a Schuler tuned vertical loop 
are complexity with the attendant increase in weight and power. To 
achieve the second integration requires typically an electro-mechanical 
integrator due to the accuracy requirements and long integration times. 

Additionally the characteristic second order servo, resulting 
from another integration in the loop, generates an oscillatory response 
when disturbed. To be useful this oscillation must be damped out. 
Proper damping requires the measurement of vehicle velocity to a 
reasonably high accuracy. 

The use of a Schuler tuned' loop is not recommended for the 
reasons stated above. Further, the increased accuracy of the vertical 
will not result in a significantly improved system since the more crit- 
ical items are the directional gyro for north direction and the accuracy 
of the distance measurement. 

4. 3 SYSTEM GAIN 

A rather arbitrary decision is made that the vertical loop 
must be fast enough, i. e. , bandwidth, such that the velocity lag re- 
sulting from a continuous input of maximum vehicle velocity shall not 
exceed 5 minutes of arc. This value is commensurate with other 
errors in the system and can be effectively tolerated. 
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Since the system proposed is a Type 1 servo we can adjust 
the velocity constant, K v> to satisfy the error at maximum angular in- 
put rate. By definition (Reference 3) 

j,, _ lim s G(s) . 

v s — ► 0 

For the system proposed (Figure' 6) 


K y = Kg 

The maximum rate of the input angle, R, is that rate cor- 
responding to a vehicle speed of 17 km/hr. 

17 

. . R = j 7 qq = .01 Rad per hour 

Therefore for a maximum tolerable error angle of 5 arc minutes 
(. 00145 Radians) we obtain 


or 


Kv = 


R 

E 


. 01 Rad per hour 
. 00145 rad 


= 6.-9 hr" 1 
Kg = 6. 9 hr" 1 


4. 4 SUITABILITY OF TIME CONSTANT SELECTION 


The value of the time constant, or loop gain, has been selected 
as 6. 9 hr” 1 . The time constant, evaluated in minutes is 8. 7 minutes 
or 520 seconds. A rule of thumb for selecting a suitable time constant 
is that the time constant of the vertical loop should be at least ten times 
longer than the longest steady state disturbance. Thus a disturbance 
lasting longer than 52 seconds would be long enough to cause a measure- 
able error in the vertical. 

The principal disturbance is acceleration, either longitudinally 
or traverse, to the vehicle direction. Consider longitudinal acceleration 
initially. It is difficult to evaluate the driver actions in relation to the 
terrain, however we should be able to get some feel for the problem 
based upon the probable distribution of obstacles that would affect the 
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pendulum. Let the vehicle stability be predominately 1/2 cps motion or 
2 sec. A pendulum with a time constant in excess of 2 sec would reduce 
the effect from vehicle dynamics and only be affected by accelerations 
that persist for longer than this period, of time. 

The length of time that these accelerations persist depends 
on the driver and the terrain. Suppose that a constant velocity is called 
for to the control system and the vehicle is to travel across the portion 
of the terrain shown in the close up photo in Ranger 7. It is difficult to 
draw a parallel pair of lines 10 feet apart, (i. e. , vehicle wheel base) 
that are not affected by the small craters. Based on this we can assume 
a maximum distance the vehicle can travel, without having to change 
direction or speed is on the order of the time represented by say^lO to 
20 feet. Even at extremely slow velocities of 1 to 2 ft/ sec this corre- 
sponds to relatively short time interval. Therefore the accelerations 
associated with these changes in velocity are rather short, on the order 
of 1 to 2 seconds. While it is theoretically possible for the vehicle to 
accelerate from a stopped position to a maximum velocity of 16 ft per 
sec at an acceleration level of . 1 ft/ sec^ (time of 160 seconds) the 
probability of its so doing is very low due to the frequent occurrence 
of obstacles. 

Another justification for the selection of the system is found 
in considering the operation of the vehicle speed control system. A 
vehicle speed control loop is a surety in the case of the unmanned auto- 
matic control phase, and likely in the manned phase of the mission. 

Thus errors in velocity are to be corrected in a short period of time, 
since the error is amplified and used to drive the motors at the maxi- 
mum rate as set by the speed control loop bandwidth. 

Further information on the vehicle and terrain characteristics 
may reveal that the low level acceleration will persist for long periods 
of time, resulting in errors in position determination. In this event it 
is proposed that acceleration measurements of the vehicle along the x v 
and y v axis be made. These can effectively be measured from body 
mounted accelerometers. The output of these accelerometers will then 
be fed to a torque r on the pendulum which will function to maintain the 
pendulum in a vertical direction during periods of acceleration and de- 
celeration. Such a technique has been used with good results (Reference 
5). 


17 


SECTION 5. 0 


LONG TERM ERROR SOURCES 


The error sources for which external compensation might be 
required are those whose frequency components are long with respect 
to the time constant of the erection loop. The following are those 
generally compensated for: 

(1) Coriolis Error 

(2) Centrifugal 

(3) Constant Velocity 

(4) Effect of rotation of the moon 

(5) Turning Errors 

(6) Component 

The magnitude of these quantities was estimated by Perkins 
(Reference 4). The Coriolis and Centrifugal errors were considered 
negligible, while correction was suggested for the velocity error, and 
latitude error. Those due to long term acceleration and turning of the 
vehicle may be ignored if the time constants for these maneuvers is 
short in comparison to the servo response time constant. 

We can now examine the magnitude of the errors due to error 
sources 3 and 4 above. 

Velocity Error Magnitude 

The transfer function for velocity error to tilt angle is derived 
from the block diagram. 

0-1.1, Ks _ 1/R 

V ” sR X K X s + Kg s + Kg 

which has a final value, 0 O = V Q /KRgas derived above. 
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The magnitude of Kg was estimated as 6. 9 hr - ^, 
a vehicle of V 0 = 5 Km/hr = 4. 56 ft/ sec. 


4. 56 ft/ sec 


(. 00192/ sec)(1728)(3381 ) ft 


= 4. 08 x 10"^ Rad 


Therefore for 


-4 


= 1.4 minutes of arc 

Thus this error need not be compensated for and can be considered in 
the error analysis as a contributing error source. 

Effect of Relation of the Moon 


The last term which may affect the accuracy of the vertical is 
that due to the latitude correction due to spin rate of the moon. The 
maximum sensed spin rate is . 536 o/hr which occurs at the lunar equator. 
The effect of this rate is equivalent to the velocity error, i.e. , a vehicle 
moving at a rate corresponding to . 536 o/hr would generate a tilt error 
of this magnitude as a maximum. 


Sine e 0 


V 

R 


V 


R(. 536 o/hr) = 


.536 o/hr (1728) (3281) ft 
57.3 o/ rad (3600 sec/hr 


14.7 ft/ sec 


- (14.7) (57) (60) 
o " (.0019*) (1728) (3281) 


4. 63 arc minutes 


It is suggested that an error source of this type be removed 
from the system by compensation. The compensation method is simple. 
The only input necessary for determining the rate of torquing the gyro 
is the latitude of the vdiicle. 
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SECTION 6.0 


ERROR MODEL' 


The coordinates of the error are taken to be along local vertical 
coordinate system axes with the x axis directed along the great circle path 
of the vehicle, the y direction along the direction transverse to the path, 
and the z axis in the direction of the local vertical. Thus a component 
of error exist along each of these axes, a track error, ex, cross track 
error, ey, and vertical error, ez. 


The first portion of the error model will be to determine the 
error angles of the reference coordinate system. The error angles re- 
presents a tilt of the reference axis in a direction given by the vector 
rotation of the axis. Thus an error angle (Cy exist due to the motion of 
the vehicle along a great circle path in the x direction. The (Dy error 
angle defines the deflection of the true vertical in the pitch plane of the 
vehicle path and results in an altitude postion error, ez. Likewise a 
error (Dz exist which relates the drift of the azimuth reference and re- 
sults in a cross track error, ey. The manner in which these error angles 
generate position errors is shown in Figure 8. The second order error 
sources which involve the product of error quantities are neglected. 

The determination of the error angles is made due to component 
error sources. Those error sources considered are as shown below 

Error angle 0/ 

Accelerometer scaling error, aj x 
Accelerometer bias, a QX 
Initial platform tilt, 0 XO 
Gyro drift rate, € y 
Error angle (Dz 

Initial alignment error, (D zo 
Gyro drift rate, ^ z 
Body spin rate, Q, 
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The error angle 0x is not determined inasmuch as its con- 
tribution is of a second order nature as shown in Figure 8. 


6. 1 VERTICAL LOOP COMPONENT ERRORS 

A block diagram of the vertical loop showing the manner that 
the component errors affect the error angle, )Dy,is shown in Figure 7. 

In appendix A, the transfer functions relating the error angle to each 
component error is derived. The results are shown in Table 2, using 
Laplace transform notation. 


TABLE 2 

RESPONSE FUNCTION FOR TYPICAL ERROR INPUTS 


Item 

Input 

— 
Transfer Function 

Input 

Response 

Vehicle 

Velocity 

V(s) 

1/R 

V 

o 

v 

°/R 

s + Kg 

s 

s (s + Kg) 

Accelerometer 

Bias 

a 

o 

K 

s(s + Kg) 

a o 

s 

a oK 

s(s + Kg) 

Gyro 

Drift 

m 

1 

MM 

£ o 

s + Kg 

s 

s(s+ Kg) 

Initial 

Tilt 

<*o 

1 

»o 

Oo 

s 

s + Kg 

s 

S ( s + Kg) 


The effect of a scale factor error in the accelerometer results 
in a gain change in the system which affects both the time constant and 
the tracking accuracy. To determine the effect of a gain change it is 
necessary to find the error angle at the correct gain which will be sub- 
tracted from the error angle with the scale factor error. This is shown 
in Appendix A. From Appendix A and Table 2 the expression for the 
total error 0y(t) may be written as the sum of the response due to each 
individual error source. 

Converting to Laplace transform to the time function the follow- 
ing expression results for 0y(t). 
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Q>y(t) = _!o_ ( !. e ~ K « t ) + 6y_, . "Kgt^ 

g ’ Kg ' 1-e ) 


+ ® 9 e ' Kg ‘ + < b * K + ^ J v 0 a 


-Wj t 


T 


a, v 
1 o 


Rw, 


- ( K + — ) v e 

v RKg ' v o e 


-Kgt 


(1) 


Where 


b = aj + 1 


= b Kg 


and the vehicle velocity is asumed constant. The last three terms re- 
present the error angle due to a scale factor error. 


6. 2 HEADING COMPONENT ERROR EFFECTS 

The system proposed utilizes a directional gyro north re- 
ference. The directional gyro may be a two degree of freedom separately 
housed instrument or a single degree of freedom gyro mounted on a stable 
element. In either case the spin axis is to be maintained in the horizontal 
plane by means of the vertical gyro. For the separately housed instrument, 
slaving signals must be supplied from the vertical gyro to the plat- 
form. It is recommended that the single degree of freedom instrument 
be utilized since the accuracy obtained from this type of instrument is 
superior to that from a two degree of freedom device. (Reference 6). 

The response function for the 0 z error angle is derived 
in reference 2. The results is equally applicable to the case studied 
here. The resulting response function for 0z(t) is given by 
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Note that for a lunar vehicle considered the term 


v 0 t 


1C 


is small (r* 1 rad at 10 hours and v Q = 17 km per hour) hence the ex- 
pression may be linearized by the small angle approximation. 


(t) Q) 0 , 


0 V 

^ox ox 


\ 


R 


t - feoz) t 


(2) 


The second term on the right hand side is negligible. Thus 
for the time and velocities considered the expression for 0 z (t) shows an 
increasing error with time. 


6. 3 Nature of the Position Error 

Assume that a vehicle is traveling along a path chosen by 
the vehicle driver. The distance that the vehicle is from the starting 
point is given by x, say, where x(t) = v(t) ’ t where 

x (t) is the distance at any time, t 

v(t) is the velocity function of time 

t is the time of computation. 


At, M 2 A-t 3 

A t 1 t— -t 1 1— 


- 

Z'H 1 1 B 


Length AB represents the actual path of the vehicle and re- 
presents the distance along the vehicle track. Let the path be broken 
up into n segments each At apart. 

Now as the vehicle is traversing the length AB the vehicles 
computer is continuously computing its position with respect to the start- 
ing point. This computation is carried out for the cross path directions, 
z and y. 

If path direction is given as along a heading \jj then the com- 
puter heading is given aslp+ 0. Then the error in the computed position 
is a function of 0, the error angle. As shown in Figure 8 there exist a 
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relationship between the error along each axis and the subsequent posit- 
ion error. This is expressed as 

ey = (0 Z ) x(t) 

ez := (0 y ) x(t) 

Now let 0 Z = 0 z (t) a function of time and x(t) = v(t) • t where v(t) is a 
function of time 

r> 

ey = Z. 0 z (t) * v(t)i • A t. 

c»*/ 1 1 

where A t are the increments of time along the path 


v(t) j the velocity during increment i 


and <Z> z (t) . is the error angle during increment i 


By definition 


£ 


<2> (t) * v(t) dt = 


lim 


h 

T 

c s > 


0 z (t) { * v(t)i -A t. 


For any arbitrary velocity profile and error angle, both expressable as 
functions of time, the total error is given by 


E y = f 0 z (t) * v(t) dt 

E z = I* 0 y (t) • v(t) dt 


The velocity of the vehicle is expected to vary dependent on the terrain 
and driver actions. However it is assumed that a vehicle speed control 
loop will be employed. Such a loop will tend to keep the velocity of the 
vehicle bounded within certain limitations, this being a matter of how 
accurately the speed control loop is maintained. It is considered a good 
assumption to let the velocity v(t) be given as 


v(t) = V Q + Vj(t) 

where v Q is the constant value of velocity about the mean value, vj (t) is 
then a function of vehicle speed and frequency of occurance of speed dis- 
trubances such as slopes. 

From photographs of the lunar surface that have been made 
it can safely be assumed that a high frequency of occurance of slope changes 
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ez = (0 y ) ( x t ) + (ey) (0 x ) 
ey = (0 z ) (x t ) + (e z ) (0 x ) 
ex = (0 z ) (e y ) + (e z ) (0 y ) 

RELATION BETWEEN ANGLE AND POSITION ERROR 



FIGURE 8 



will occur, say at least one every 10 to 20 meters. At a mean velocity 
of 5Km/hr (1. 39 m/sec), the slope changes causing a load on the velocity 
servo will occur- at a frequency of 

f- 1-39 m/ sec _ 139 slope variation/ sec 

10 m/ slope variation 

This is to be compared with the bandwidth of the servo loop 
which was chosen as 

6 . 9 Rad/hr = ra ? 6 . 28 cycle_(hr^_ 3 . 012 les/sec . 
lv rad 3600 sec 

Thus the perturbation for the assumed speed and distance between slope 
variations are outside the bandwidth of the velocity servo and hence to 
do affect the lag angle of the gyro. 

From another standpoint the distance traveled over a period 
of time will be equal to 

d = V Q (t) + Vj (t) 

If vj is assumed to have a statistical distribution about the steady state 
velocity v Q then the distances represented by vj ( t ; are equally probable 
for substracting from or adding to the distance computed v Q . t .*. for 
practical purposes Vj(t) may be neglected. 

For a servo controlled speed loop then the v(t) may be re- 
presented as v(t) = v Q t* 
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SECTION 7. 0 


ERROR EVALUATION 


The system evaluated is a Type 1 single degree of freedom 
vertical loop with a separate directional gyro for an azimuth reference. 
The value of the error terms and parameter values of the error terms 
is given in Table 3 below. 


TABLE 3 

PARAMETER VALUES FOR SYSTEM CONSIDERED 


It em 

1 1 

Values 

Units 

Accelerometer Bias 

a o 

10" 3 

g 

Gyro Drift 

^ y z 

.05 

degrees per hour 

. 1 

degrees per hour 

Loop Gain 

Kg 

6.9 

hr - 1 

Body Radius 

R 

1728 

kilometers 

Vehicle Velocity 

V 

5 

kilometer per hour 

o 

10 

kilometer per hour 

Scale Factor Error 

a l 

5 

percent 

Initial Tilt (vertical) 

0 

oy 

6 

minutes of arc 

Initial Alignment 
(azimuth) 

0 

oz 

2 

minutes of arc 

time 

t 

6 

hours maximum 


7.1 INITIAL ALIGNMENT ERRORS 

The initial tilt of the vertical, 0 q , is dependent on the align- 
ment method and the time allotted for the imtial tilt to converge to a 
minimum value. Since the vertical loop is equipped with a servo to erect 
the spin axis then a matter of 3 time constants (say 20 minutes) is required 
to have the loop converge to its minimum value. For the error analysis 
an initial error of 6 minutes of arc is assumed. This would correspond 
to some manual method of rotating the platform to a nominal level con- 
dition prior to energizing the erection loop. 
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The initial heading alignment is much more critical. Since 
the heading gyro operates in a free condition the initial misalignment 
will result in position errors that become significant with time and distance. 

Three basic techniques may be used to initially align the 
axis of the direction memory. These are; 

(1) Surveying (landmarks or celestially) 

(2) Gyro compassing 

(3) Radio Celestial in cooperation with CSM 

For the case of surveying by celestial observations. The 
accuracy with which an azimuth may be determined by star sightings is 
dependent on the following items: 

(a) Unknown libration of moon 1/2 min. 

(b) Accuracy of knowledge of the location of observation (1 min) 

(c) Accuracy of the reading (1/5 min) 

(e) Sight reduction (2 min) 

Considering each of these error sources as random variables 
the rss of the resultant accuracy is 

=y 1/4 + 1 + 1/25 +4 ' 

= 2.4 minutes of arc 

The second method is that of obtaining north by a north seek- 
ing gyro. Estimates have been made on the accuracy of a two axis 
spherical north seeker as 2 minutes of arc after a period of 2 hours at the 
lunar equator. 

The third method is subject to an interface problem with 
availability of the CSM. Additional equipment would also be required 
to compute the direction to the orbit ground track. Hence this third 
method is disregarded. 

In summary the initial alignment accuracy of the azimuth 
direction memory can be taken as a one sigma deviation of 2 minutes 
of arc. 

7 . 2 EVALUATION OF NOISE ERROR SOURCES 

The following error sources are considered: 
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Accelerometer Bias 


Accelerometer Scale factor 
Gyro drift rate 
Initial Tilt 

All the above are random in nature and hence the root sum of 
the squares of the error from each source will represent the desired 
accuracy. Expressed as a sum, the expression for the tilt of the ver- 
tical Oy(t) is given by equation { in Section 6-1) 

7.2.1 Vertical Tilt due to Accelerometer Bias 

The error due to accelerometer bias is expressed as: 

0 y (t)= ( l-e _at ) 

2 The results of this expression for an accelerometa: bias 
of 10 g is given in Table 4. The error is given versus time. The 
third column gives the tilt error in milliradians as one standard deviation 
{ 1 (T value.) The 1CT value is given in column 4. 

TABLE 4 


VERTICAL TILT ERROR DUE TO ACCELEROMETER BIAS 



The error angle due to gyro drift rate is expressed as 

= I 1_e " at * 
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The results of this expression for gyro drift rates of .05 
and . 1 degree per hour gyros is given in Table 5. The angle is given 
in terms of 1 <S and 1 <f 2 values for each of the drift rates considered: 


TABLE 5 

VERTICAL TILT DUE TO GYRO DRIFT RATE 


time 

(hours) 

£ o = . 05 deg/hr 

6 q = • 1 deg/hr 

TO’ 

mr 

icr 2 

KT 

mr 

icr 2 

0 

0 

0 

0 

0 

. 1 

. 126 

.0159 

.063 

.0044 

.2 

. 189 

.0357 

.0945 

.00893 

.4 

.236 

.0558 

.118 

.0139 

.6 

.248 

.0613 

. 124 

.0154 

.8 

.252 

.0635 

. 126 

.0159 

1 

.253 

.064 

. 1265 

.016 

2 

.253 

.064 

. 1265 

.016 

4 

.253 

.064 

. 1265 

.016 

6 

.253 

.064 

. 1265 

.016 


7.2.3 Vertical Tilt Due to Initial Alignment Error 


The error angle due to an initial alignment error is given by 

Oyf') = «■„ 

The values of 0 as a function of time are computed and 
listed in Table 6. ^ 


TABLE 6 

VERTICAL ERROR DUE TO INITIAL ALIGNMENT ERROR 


time 
( hours) 

r 

-at 

e 


1<T 2 

0 

i 

1.74 

3.04 

. 1 

. 502 

.376 

.778 

. 2 

.2502 

.437 

. 191 

.4 

.0635 

.108 

. 0157 

.6 

.016 

.0279 

0 

.8 

.004 

.007 

0 

1 

0 

0 

0 

2 

0 

0 

0 

4 

0 

0 

0 

6 

0 

0 

0 
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7.2.4 


Vertical Tilt Due to Accelerometer Scale Factor Error 


The scale factor error is determined for a change of 5% in 
the system gain. This error may result from environmental changes 
of the equipment and is to be expected. In this system the scale factor 
change results in some random change of the vertical and contributes 
indirectly to a position error. This error source may be thought of as 
any random dlstrubance which may cause a change in vertical loop gain. 

The error was computed for a velocity of 5 km per hour and 10 km per 
hour and was found to be completely negligible. 

7.3 EVALUATION OF 0 V SYSTEMATIC AND BIAS ERROR SOURCES 

Two bias errors sources are considered which result in tilt 
of the vertical. These are the velocity lag and the latitude error. 

7.3.1 Velocity Error 

,-i ns 1 \ -k t , V o 

^(t) ' (K ‘ Rkg ) v o e g Rkg 

The assumed velocities were 5 & 10 Km/hr. In each case 
the final steady state error will exist as o which is a maximum 
peak. KgR 

_ c . , w = 1.44 minutes = . 42 mr 

For v =5 km per hour Q. 
o t 

_ , =2.88 minutes =.82 mr 

For v =10 km per hour Q 
o t 


7.3.2 


Latitude Error 


Spin Rate of Moon = 


14.7 ft /sec km 
73381' ft) 

15. 6 km/hr 


36 00 sec 
hr 


) 


0 =4.5 minutes max = (1.28 mr ) 


7.4 ERROR SUMMATION FOR TERMS 

y 

Table 7 summarizes the errors in the following manner. 
The total peak error is expressed as 


/ r 2 2 

Total Peak Error = / 3 (Noise Error) 4 (Bias Errors) 

The results of this tabulation are shown in Figure 9. 
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SUMMATION OF ERRORS IN INDICATED VERTICAL 





















































































































7.5 EVALUATION OF 0 VERSUS TIME 

z 

The equation for 0 (t) was shown in section 6. 2. as equation 
(2). There are three principal error sources that reflect in an azimuth 
error angle. These are gyro drift, uncompensated latitude rate due to 
the spin of the moon, and the initial alignment. These are considered 
as random variables. 

7.5.1 Initial Alignment Error 

The initial alignment error of the azimuth reference 0 
is a constant value. It was estimated in section 7.1 as a one sigma value 
of 2 minutes of arc. 

0 =2 minutes of arc = .581 mr 
oz 

1 (T ^ value of 0 = . 338 

oz 

7.5.2 Azimuth Gyro Drift 

The azimuth gyro drift is estimated at between .05 and . 1 
degree per hour. The assumption is made that these drift rates are 
constant. The resulting error angle after a period of time for the stated 
drift rates is shown in Table 8. 


TABLE 8 

DRIFT ANGLE VERSUS TIME FOR AXIMUTH GYRO 

2 

1 Values, 0 Milliradiaus 1 Values 


1 

(«<) <*) 


0 milliradiaus 
z 

time 

(hours) 

v = 5 km/hr 

o 

v = 10 km/br 
o 

v =5 km/hr 
o • 

v q = lOkm/hr 

0 

0 

0 

0 

0 

. .1 

.0872 

. 1744 

.00762 

.0304 

.2 

.1744 

.3488 

.0304 

. 122 

.4 

.3488 

. 6976 

. 122 

.487 

.6 

. 5232 

1.0464 

.273 

1.09 

.8 

.6976 

1.3952 

.487 

1.94 

1 

.872 

1.744 

.761 

3.04 

2 

1.744 

3.488 

3.04 

12.22 

4 

3.488 

6.976 

12.22 

48.7 

6 

5.232 

10.464 

27.3 

109. 
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7.6 


ERROR SUMMATION FOR 0 TERMS 

' 11 Z 

Table 9 summaries the two contributing errors. The total 
peak error in the indicated azimuth is three times the v. s. s. sum of 
the two independent random errors. The results of the tabulation is 
shown in Figure 10. 
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SUMMATION OF ERRORS IN INDICATED AZIMUTH 

































































FIGURE 10. TOTAL PEAK ERROR OF INDICATED AZIMUTH REFERENCE 
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7.7 


DETERMINATION OF ALTITUDE POSITION ERROR, ez, 
FROM NOISE SOURCES 


The position error is to be calculated for each error source 
in turn, utilizing the formula given in section 6. 0. The total position 
error for the noise sources is then given by the r. s. s. sum of the error 
due to each individual error source. The error contribution for each 
source is given in the following paragraphs. 

7.7.1 Position Error Due to Accelerometer Bias 

The tilt angle for an accelerometer bias error, a , for the 
assumed system, is given by. l " 

0 y (t) = ^ ( 1 -e " Kgt ) 


The position error is given by 



The values of e versus time are shown in Table 10 and presented 
graphically in Figure 11. 

TABLE 10 

POSITION ERROR DUE TO ACCELEROMETER BIAS 


time 1 ~ 

(hours) 1 + Tci" e " Kgt V o = 5 km/hr v c = 10 Whr 



meters 

meters 

. 145 

0 

0 

. 1728 

. 139 


278 

.2363 

.456 


913 

.4092 

1.3210 

2. 

64 

.602 

2.285 

4. 

57 

.800 

3.275 

6. 

55 

1 

4.275 

8. 

55 

2 

9.275 

18. 

55 

4 

19.275 

38. 

55 

6 

29.275 

58. 

55 
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FIGURE 11. ALTITUDE POSITION ERROR DUE TO ACCELEROMETER 
BIAS (Constant Velocity assumed) 
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7.7.2 


Position Error Due to Gyro Drift 


The tilt angle for the vertical gyro, due to gyro drift, 
for the assumed system is given by: 

V‘> ■ Sr < 1 - ' Kgt > 

The position erro_r is given by 

e z = a / V o V t>,dt 

which reduces to 

v o €p I" , 1 ^ -Kgtl _ V Q € o 

6 z kg Kg J (Kg) 2 

Tie values of e versus time are shown in Table 11a and lib plotted in 
Figure 12. 


TABLE 11a 


POSITION ERROR DUE TO GYRO DRIFT ( € Q = . 1 o/hr) 



€z * 

meter 

time 
(hour s ) 

t 1 1 c " Kgt 

c 

o 

II 

Ui 

v =10 
o 

* Kg e 

0 

.145 

0 

0 

.1 

. 1728 

.035 


.2 

.2363 

. 115 

.232 

.4 

.4092 

. 335 

. 670 

. 6 

.602 

. 578 

1.156 

.8 

.800 

.826 

1.652 

1 

1 

1.084 

2.168 

2 

2 

2.446 

5.892 

4 

4 

4.876 

9.752 

6 

6 

7.396 

14.792 
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FIGURE 12. ALTITUDE ERROR DUE TO GYRO DRIFT 
(Constant Velocity Assumed) 



TABLE 11b 


POSITION ERROR DUE TO GYRO DRIFT ( = .05 o/hr) 


time 

(hours) 

. , A :_._. C -Kgt 

£ Zj , Meters 

v =5 
o 

v = 10 
o 

* + K g 

0 

. 145 

0 

0 

.1 

. 1728 

.017 

. 135 

.2 

.2363 

.057 

.115 

.4 

.4092 

.1675 

.335 

.6 

.602 

.289 

. 578 

.8 

.8005 

.413 

.826 

1 

1 

.542 

1.084 

2 

2 

1.223 

2.446 

4 

4 

2.488 

4.876 

6 

6 

3.798 

7.396 


7.7.3 Altitude Position Error Due to Initial Tilt 


The initial alignment error of the vertical gyro will converge 
when the loop is energized. The effect of the initial tilt error will cause 
a position error that is constant with time. The expression for the 
error angle and subsequent position error are as given below 


V‘> ■ *o * ‘ Kgt 


e 

z 




The position error resulting from this error source is shown in Table 
13 & Figure 13. 
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TABLE 12 


POSITION ERROR DUE TO INITIAL ALIGNMENT TILT 


time 

(Hours) 

1 e - K & 
Kg 

E , Meters 
z 

v =5 km/h 
o 

r v — 10 km) /hr 

o 

0 

. 145 

0 

0 

.1 

.0728 

.91 

1.82 

.2 

.0363 

1.371 

2.742 

.4 

.0092 

1.613 

3.326 

.6 

.00232 

1.8 

3. 604 

.8 

.00058 

1.82 

3. 64 

1 

0 

1.83 

3. 66 

2 

0 

1.83 

3. 66 

4 

0 

1.83 

3.66 

6 

0 

1.83 

3 . 66 


7.7.4 Position Error Due to Scale Factor Error 


The position error resulting from a scale factor change of 
five percent was found to introduce no significant position error. 

7.7.4 Position Error Due to Readout Error 


The position computer will calculate the vehicle's position 
with reference to the angle on the readout. It is assumed for this sytem 
that the readout accuracy is correct to within 3 minutes of arc. This 
will also be a random function. 


The error in the readout angle is denoted-0 R . has a 1C- 

value of 3 minutes of arc = .873 milliradians . The position error e 
is given by z 


e 

z 
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FIGURE 13. ALTITUDE ERROR DUE TO INITIAL TILT ERROR 
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TABLE 13 


POSITION ERROR DUE TO READOUT ERROR 



7.7.5 Position Error Due to Vehicle Acceleration 


The effect of vehicle acceleration and deceleration is to deflect 
the pendulum. This will result in an error in the indicated vertical which 
affects a position error. The magnitude of this postion error is determined 
in Appendix B for typical acceleration levels. The position error is 
generated during the portion of time that the acceleration level persist, 
and the total error becomes the summation of the errors from each 
acceleration or deceleration period. 

It is recommended that the acceleration generated errors be 
compensated for. It is a simple matter to measure vehicle acceleration. 
Tlis acceleration is then to be used to torque the gyro back to the true 
vertical, as done in the system in Reference 5. 

7.8 ALTITUDE POSITION ERROR RESULTING FROM BIAS ERRORS 

Altitude postion errors due to bias sources arise from the 
velocity lag for a type 1 vertical loop, and due to the Corolios error. 

The Corolios error is to be removed by compensation. 

Velocity Error 

The expression for the error angle due to an initial step velocity 

input is 

0 (t) = ^g— ( 1-e _Kgt ) 
g 


46 








The resultant expression for the position error becomes 


e 

z 


v 


o 

= RKg 


(t + 


1 

K g 



o 

R (Kg) L 


The value of this expression for the assumed case is given in tabular 
form in Table 14 and graphically in Figure 14. 


TABLE 14 

POSITION ERROR RESULTING FROM VELOCITY LAG 




E Error, Meter 
z 

Time 

(Hours) 

<‘*K g ^ 

<J 

o 

li 

v = 10 
o 

0 

. 145 

0 

0 

. 1 

. 1728 

.059 

.232 

.2 

.2363 

. 192 

.764 

.4 

.4092 

. 556 

2.22 

.6 

.602 

.961 

3.84 

.8 

.800 

1.376 

5. 50 

1 

1 

1.796 

7.18 

2 

2 

3.896 

15. 58 

4 

4 

8.096 

32.38 

1 

6 

12.296 

49.18 


Gravity Bias Error 

The deflection of the vertical due to gravity anamolies may be 
considered as a bias error rather than a noise error because of the 
limited range of operation of the vehicle in relation to the moon's size. 
The period of operation of the dead reckoning function is approximately 
six hours or less and covering an area not greater than 40 km across. 
Over this area, large measureable vertical deflections, as may be cause< 
by mountain masses, will likely remain constant. 

The total vertical deflection will be distributed along the 0 
direction or the 0 direction dependent on the heading of the vehicle. 

The maximum one^ sigma deviation of gravity caused vertical errors has 
been assumed as one arc minute. Ijf the heading of the vehicle can be 
assumed as an arbritrary random function then the deflection along the 
0 axis, or the 0 axis, can be assumed as the rms values of the sine 
fiSiction, or in tins case, as .707 arc minutes. 
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FIGURE 14. ALTITUDE POSITION ERROR DUE TO VELOCITY LAG 


48 




















The position error resulting from this source, governed by 
the assumptions made above, is given by 

e = v (t) (0 ) 

z o Yg 

where 0 is the deflection of the vertical due to gravity effects. Table 
15 lists this error as function of time 


TABLE 1 5 

ALTITUDE ERROR DUE TO GRAVITY AN AMO LIES 



e Error, Meters 
z 

Time 

(Hours) 

r km 

v =5 

o mr 

v - 10 ^ 

o hr 

0 

0 

0 

.1 

. 103 

.206 

.2 

.206 

.412 

.4 

.412 

.824 

.6 

.618 

1.236 

.8 

.823 

1.648 

1 

1.03 

2.06 

2 

2.06 

4. 12 

4 

4. 12 

8.24 

6 

6. 18 

12.36 


7.9 ERROR SUMMATION FOR ALTITUDE ERROR 


Table 16 and 17 contain the results of an error summation of 
the various error sources into an overall peak error for an assumed 
velocity of 5 km/hr, and 10 km/hr. The results of this are shown 
in Figure 1 5. 
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TABULATION FOR ALTITUDE ERROR, E , v = 5km PER HOUR 
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SECTION 8. 0 


DETERMINATION OF CROSS TRACK POSITION ERROR, ey 


The cross track error is computed in the same manner as the 
altitude error. The position error is given by 



8. 1 CROSS TRACK ERROR, ey, DUE TO INITIAL ALIGNMENT 

ERROR 


The cross track error is given by 

0z = 0 O z 

The resulting position error becomes 
0oz v Q t 

This error is tabulated in Table 18. 


TABLE 18 

CROSS TRACK ERROR, ey, DUE TO INITIAL ALIGNMENT ERROR 

Error ey, Meters 


Time 

(Hrs) 

v 0 = 5 km/ hr 

v Q = 10 km/hr 

. 0 

0 

0 

. 1 

. 291 

. 582 

, 2 

. 582 

1. 165 

. 4 

1. 165 

2. 33 

. 6 

1. 75 

3. 5 

. 8 

2. 32 

4. 66 

1 

2. 91 

5. 82 

2 

5. 82 

11. 65 

4 

11. 65 

23. 3 

6 

17. 5 

35 
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8 . 2 


CROSS TRACK ERROR, ey, DUE TO GYRO DRIFT 


The cross track error angle is given by 
= ^oz * 


The resulting position error is 

re 


ey(€) = 


v o foz t dt 


Vo 


e G t c 


The value of this error is given in Table 19. 

TABLE 19 

CROSS TRACK ERROR, ey, DUE TO GYRO DRIFT, fc QZ = . 05°/hr 


Time 

(Hrs) 

ey error, meters 

vo - 5 km/ hr 

v Q = 10 km/ hr 

0 

0 

0 

. 1 

. 0436 

. 0872 

. 2 

. 1745 

. 348 

. 4 

. 7 

1. 395 

. 6 

1. 58 

3. 14 

. 8 

2. 79 

5. 58 

1 

4. 36 

8. 72 

2 

17. 4 

34. 9 

4 

69. 7 

139. 5 

6 

157. 

316. 


8. 3 CROSS TRACK ERROR, ey, DUE TO READOUT ERROR 

The position computer will calculate the vehicle's position with 
reference to the angle on the azimuth readout. It is assumed for this 
system that the readout accuracy is correct within three minutes of arc. 
This will also be a random function. 

The error in the azimuth angle readout is denoted 9^. 9 a has 

a one sigma value of three minutes of arc = . 873 milliradians. The 
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position error is given by 


f Q A v o dt = e A v Q t 


This error is listed in Table 20. 


TABLE 20 

POSITION ERROR, ey, DUE TO READOUT ERROR 



ey Error, Meters 

v 0 = 5 km/hr 

v Q = 10 km/hr 

0 

0 

. 436 

. 873 

. 873 

1. 745 

1. 74 

3. 59 

2. 62 

5, 24 

3. 49 

6. 98 

4. 36 

8. 73 

8. 73 

17.45 

17. 4 

35. 9 

26. 2 

52. 4 


POSITION ERROR SUMMATION FOR CROSS TRACK, ey 


ERROR TERMS 

Tables 21 and 22 summarizes the major contributing error to 
cause a cross track error. The total peak error is three times the 
r. s. s. of the individual error sources. The resultant error is plotted 
in Figure 16. 

TABLE 21 

SUMMATION OF ERROR TO AFFECT CROSS TRACK ERROR, 



0 . 0848 . 339 1. 36 3. 06 5. 39 

0 ^ 0 .029 .49 2. 5 7. 8 

0.19 . 77 3. 04 6. 88 12. 2 

0 . 9 08 1. 85 3. 83 6. 11 8. 72 



34 136 306 

304 4870 24600 

762 306 688 

34.12 126.24 277 
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TABLE 22 


SUMMATION OF ERRORS TO AFFECT CROSS TRACK ERROR, 

v Q = 10 km/hr 






i <r 2 

Values 










Time - 

Hour 







0 

■ 

. 2 

• 4 

. 6 

. 8 

1 

2 

4 

6 

Initial 

Alignment 

0 

.339 

1. 36 

5.41 

13. 25 

21. 7 

33. 9 

136 

543 

1225 

Gyro Drift 

0 

A J 0 

. 121 

1. 96 

9. 9 

31. 1 

76 

1220 

19400 

100000 

Readout 

0 

.77 

3. 04 

12. 8 

27. 5 

49 

76. 2 

306 

1290 

2750 

Total Peak 
Error 

0 

1 

1. 82 

3. 69 

7. 77 

12. 32 

17.48 

1 

23. 48 

70. 61 

252.4 

558. 5 
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SECTION 9. 0 


CONCLUSIONS AND RECOMMENDATIONS 
9. 1 CONCLUSIONS 

A navigation system of the type discussed herein will provide 
sufficiently accurate position information to allow safe conduct of a 
lunar vehicle for missions involving times up to six hours. Beyond six 
hours, the error in the track and cross track directions for typical 
velocities of 5 to 10 kilometers per hour will become excessive; on the 
order of 300 meters in the track direction (due to distance measure- 
ment inaccuracy) and 300 meters in the cross track direction (due to com- 
ponent inaccuracies). This error is sufficiently large that the error 
in absolute position of the vehicle approaches one kilometer. This is 
obtained when combining the dead reckoning error with the error in the 
initial position fix. At this point a new position fix should be made. 

The theme of this investigation was to study a relatively simple 
system configuration, not relying on high accuracy gyros and complicated 
compensation methods. The drift rate of the azimuth gyro was selected 
as . 05 degrees per hour. It is not anticipated that a small, lightweight, 
rugged gyro, with higher accuracy, suitable for use on the moon will be 
available during the time frame of the first lunar exploration. 

Relatively simple vertical gyro instrumentation will achieve 
altitude computation to an accuracy of 70 meters. 

The acceleration compensation is the singlemost critical item 
affecting system error. If the acceleration error induced by the vehicle 
can satisfactorily be compensated for then a system of the sort discus- 
sed in this report is entirely satisfactory. If the acceleration com- 
pensation cannot be mechanized then a more complicated system, i. e. , 
Type II vertical loop leaves the only alternative for the maintaining of 
an accurate vertical. 

9. 2 RECOMMENDATIONS 


It is recommended that a dead reckoning system be bread- 
boarded and tests performed on its accuracy. This should be in cross 
country terrain situation to simulate vehicle accelerations that may be 
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encountered on the moon. This will allow the determination of the 
velocity and acceleration profile of a vehicle. When this information 
is available an analysis may be performed on a computer to determine 
the best system. 
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APPENDIX A 


RESPONSE FUNCTION OF VERTICAL LOOP TO TYPICAL 

INPUT ERROR SOURCES 


From the block diagram, Figure 8, the loop gain function 
for each input is obtained. 


(a) 


(b) 


Velocity Input, v q 

Qyjs) _ l _ s 

V ( s ) 1 + Kg s + Kg 

RS s 

V(s) = v q /s 
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(d) 


Cj)y(s) = 




s + Kg 
Scaling Error, a^ 


The error angle in the presence of scaling error is deter- 
mined as follows: 
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b 2 k 

Oy(s) = s V(s) y + b 2 Tg- 


V(s) 
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s + w. 


s + w, 


letting b = + 1 


^ , v _ s V(s) ai K + s Y(s)K - V(s) /R 

y y (s) = — — 


s + w 


1 


The error angle with the correct gain is 


a i-s - s V(s) K - V(s)/R 
Oy(s) rTKg 

Thus the error angle due to scale factor error becomes 

_ s V(s) ai K + s V(s) K - V (s)/R _ s V(s) K - V(s)/R 

^(s) s + Wj s + Kg 


(e) 


Gyro drift , 
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APPENDIX B 


RESPONSE FUNCTION OF TYPE I VERTICAL GYRO SYSTEM 
TO A CONSTANT ACCELERATION 

The error angle response of a type 1 vertical loop is given by 

«•> = K ° ; |/ e R v(b) in 

a 

Let V(s) = ^ ■ for a constant acceleration, equivalent to a ramp input 

s 

of vehicle velocity. 


0(s) 


Ks - 1/R 
s + Kg 


a 

o 



0(s) 


a v 
o K 

s(s + Kg) 


a o/R 

s 2 ( s + Kg) 


( 2 ) 


Evaluation of 1st term in equation (2) by partial fraction expansion: 
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B = a-* o (s+ Kg) 
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(4) 
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Evaluation of 2nd term in equation (2) by partial fraction expansion. 


a o/R 

s 2 (s + Kg) 



C 

~+ Kg 


(7) 


b3 


multiply by s 


a o/R _ A + Bs + Cs 2 
s + Kg “ s + Kg 


(8) 


let 8=0 



(9) 


differentiate equation (8) with respect to 8 


a o/R 0 _ B + 2Cs 
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(s + Kg) 


( 10 ) 


let s - O B = R ( K g ) £ 


(ID 


multiply equation (7) by (s + Kg) 
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2 2 s 
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let s = - Kg 
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a 

o 

IT 
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substituting equations (6) and (14) into equation (2) the Laplace expression 
for the total 0 (s) is obtained. This then reduces to the following time 
expression for 0 (t) due to a constant acceleration. 
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(15) 


®<«> = f- ( i - . - Kg j 


+ <*> + rJKI )^ 1 ' e ’ Kgt ) 


RKg 


The first and third terms of equation (15) reduce to constant values after 
several time constants, leaving the dominate error as term two which 
increases linearily with time. 

Position error resulting from a constant acceleration . 

The position error is given by 

e z = Jt>m v(t) dt (16) 


0 (t) is given by expression (15) 


v(t) for a constant acceleration = a t 

o 


• • e z 


[2 2 

1 a a 

o ■ L o 

— l^gT 2- 



*) ( 1-3 " kgt ) dt + 0 


RKg 



dt 

(17) 


r 0 JO 

Evaluation of this expression between the limits of o and t results in 
the following expression for the position error. 


e = a 
z o 




g + RTKg) 2 J\^7 + 3RKg y + TRKg ( (18) 

^ ^ J 2 ; 

Tlis term is evaluated for three acceleration levels, .5 ft/sec , 1 ft/sec 

and 2 ft/sec . The maximum velocity capability is assumed to be 20 ft/ s 

Therefore, maximum times for evaluation are, respectively, 40 sec, 20 s 

and 10 sec. The results of this evaluation is shown in Table 1 . 
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TABLE 1 


ALTITUDE POSITION ERROR VERSUS TIME 



(Type 1 vertical loop with loop gain, 
for constant acceleration) 

Kg, = 6. 9 hr' 1 , 

Time 
(Sec . ) 

Resultant Postion Error 

a = . 5 ft per sec^ a 

o o 

in feet 
= 1 ft per 

sec^ a Q = 2 ft per sec 

5 

151. 5 

606 

1212 

10 

306 

1225 

2450 

15 

457 

1828 

3656 

20 

606 

2425 


25 

756 



30 

903 



35 

1060 



40 

1196 
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